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ABSTRACT 

The practice of utilizing mathematics instruction, in 
which only^practical engineering applications and not the pure 
mathematical proofs are stressed, is questioned • Three cases Sire 
presented in which students made erroneous conclusions concerning an 
engineering topic. In each situation, the student violated a 
mathematical principle which was not evident to him because his 
instruction had dwelled only on the practical applications 'of the 
formulae. , (CP) 
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MATHEIviATICAL FORiMALITIES AW EFGirEERING ^STCDIES . 

Many engineering schools have incorporated a mathematics 
department in their units with the explicit objective to 

teach "engineering mathematics", as opposed to pure or rigor- 

I 

ous mathematics. The subjects in mathematics are taught for 
direct application and "unnecessary or tedious" proofs are 
omitted or at best presented a^" unimportant . As a conse- 
quence, engineering students do not have much opportunity to 
acquaint themselves .with carefully presented mathematical de- 

veloTDments . 

* *■ • 

Since the engineering applications presented in class 
are already well developed and tacitly assumed to be cor- 
ect, it may appe^ar that no need exists to proceed carefully 
through mathematical proofs. The instructor omits such 
proofs in his lectures, often unintentionally "to save.ti^e", 
or^sometimes intentionally because he considers them "irrele- 
vant" mathematical formalities. For example, no explicit 
requirement may be state^^that before dividing by a factor, 
this factor should nol; be' zero, 8r no mention may be madai :r 
that a series expansion may be convergent onl^ in a certain.,^ 
range . 

-,The objective of this pap'er is to discuss personal 
teaching exDeriences in which the lack of a basic under- 
standing of mathematical principles precluded a clear 



comprehension of some engineering concepts, ^ 
The consequence of too much confidence ! 

r 

Too much confidence in "proven" engineering methods may • 
lead to inadmissible errors as demonstrated in this, first 
case presented. In the finite element method for the dynamic 
analysis of framed structures, the axial, torsional, and 
flexural actions are expressed in terms of 4;he elastic and 
inertial forces constituting respectively the stiffness and 
mass matrices. These matrices relate* the forces to the < 
displacements or to the accelerations of a basic beam element 
and are derived from assumed disp;^acement functions. Specif i- 
-cally, it is assumed that the static deflection curve of a 
beam element also represents the deflection for a^ dynamic 
situation* This approach permits discretizing of the continuous 
member by. means of its elastic properties as given in the 
stiffness matrix, and its inertial characteristics as given 
in the consistent mass matrix^. A more crude and simple method 
of disjJretizing the inertial properties of a beam, element is 
to simply allocate half of its total mass to each end as a . 
concentrated lump mass.' The discretizing process has also 
been carried o\x% mathematically by means of the apprqpiate 
differential equation. From this is derived the dynamic stiff- 
ness matrix -which expresses the exact relation between harmonic 
forces and displacements in terms of transcendental trigonometric 
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. and hyperbolic functions. 

Some time ago I proposed as a topic -for a Master's Thesis 
a parametric comparative study of the \hree methods of discre- 
tizing a beam element. The student working on this subjiect 
wrote a computer program and proceeded to test the program 
selecting simple numerical values. Specifically, he gave unit 
values to the modulus of. elasticity, moment of inertia, a^ 
otlier parameters. It happened that all these three methods 
which normally r>roduce relatively close results, for the 

- r 

simple unit values gave widely different answers. This 

unexpected result certainly added more interest to the subject 

2 

of this thesis and resulted in a technical paper in which it 
v;as shown ^ that so-called stiffness and mass maTtrices v/ere ' 
nothing other than the first two terms of the pov/er series 
expansion of the exact dynamic stiffness matrix. The mathe- 
matical approach Used to obtain these matrices, in contra- 
distinction to the' engineering approach, allowed the specifi- 
cation of/ the range convergence of the resulting series. It 
was found that. the unitary^numerical values used in testing 
the computer program were outside the range of convergence of 
the series. Consequently, the fact that these values led to 
divergent results with large errors was no longer surprising 
or unexDected. 
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Make the. determinant equal zero ! 

The experience described above led to the following second 
case. As has already been stated the use of the exact equation 
fo^ the discretization of framed structures yields the dynamic 
stiffness matrix containing trigonometric and hyperbolic func- 
tions. The common procedure to find the natural frequencies 
of such structures is to set the determinant of the dynamic 
stiffness matrix equal to zero. In the search for. zeros, the 
determinant is evaluated for a series of values of frequency 
in the region of interest. A change in the sign of the value 
of the deteiTninant indicates the presence of one ..or more natural 
frequencies. In other words, to find the natural frequencies, 
the analyst sets the determinant equal^ to zero and searches 
for its roots. As it happens, the searching process is not 
possib^ for certain values of the argument; ' far these values 
the determinantal function is not defined. Thi%^ difficulty • 
occurs in the study of flexural vibration as v/ell as for tor- 
sional or axial vibration. The situation can be ^illustrated - 
for the simple case of an axially-loaded rod.. The governing 
equation in this case, as given in referenc-e'^, is the v/ave 
equation ' . 

a^u - _m = 0 . (1) 

8 AE 8"t2' ' - 
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.where m is the mass per unit of length, A the cross-^sectional 
area, E the modulus of elasticity, and u the displacement at 
coordinate x and time t. 

The solution of equation (l) yields a harmonic motion of 
amplitude 

u = sin 0X + cos gx - " (2) 

where 



(Jt^is the natural frequency and C^, Cg^are constants , of inte- 
gration. ... , 

To obtain the dynamic stiffness matrix for the axially ^ 
vibrating' element, boundary conditions indicated by equations 
(4) and '(5) are imposed. 

u{x = 0) 6^ ■ u{x = L) = (^) 
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o 

du(x = 0) = - PA^ . du(x = L) = PB^ 

dx . ' • • AE ' • . dx - AE 

where 5^, 63 and Pg are respectj.vely the displacements 

and fflrces in the beam element as shown in Fig. 1.'- 

+ A,E . ^ ■• 

U- - L ; 'J 

Pif:. 1- Displacement , and forces for axlally loaded beam 
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Substituting the boundary conditions of equations (4) and (5) 
into equation (2) results in 



and 



~ 6 " 
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singL 


cos3L 


■"Pa "' 




-1 


0 




= AEB 








Pb . 




COS3L 


--sinBL 



'2 _ 



C2 



(6) 



(7) 



Equation (6) is solved for the constant of integration, 
namely 



J 



-cotBL . GOsecgL- 
1 0 



-«B 



(8) 



subjected to the condition 

sin&L 0 



(9) 



Subtitution of eq,uation (8) into equation (?) fives equation 

(lO) which relates -the forces to the dispJ^c^ments at the 

* . ■ 

coordinates through. the dynamic matrix. 



-Pb J 



EAB 



cotgL -cosecgL 
-cosecBL cotBL 





^«a' 







(10) 



It is interesting to observe that it .is. not possible to 
find the natural frequencies for an axially loaded beam by 
simply letting the determinant of the dynamic stiffness matrix 
in equation (10) be equal to zero; Subjected to the condition 
given by equation (9)i the determinant of this matrix has a* 



constant .value, minus one. Furthermore, the condition ^^iven 

by equation (9)-» v;hich is required to obtain equation (10), 

precludes the determination of * the: Jiatifral. frequencies from 

the determinantal function. The zeros of sin are precisely 

the natural frequencies of a beam element undergoing axial 

« 

vibrations. For such an element, the natural frequencies may 
be found by equating to zero the* determinant of the 2x2 
matrix in equation (7)t but naj: by using the dynamic stiffness 
matrix as has been demonstrated. This difficulty arises 
because the natural frequencies of the system formed by only 
one member or^ element are obviously the same as the natural 
frequencies of that member. In an actual problem where^'fKe^ 
dynamic stiffness matrix for the system is obtained from 
several elements, the determinant of this matri-x will not be 
defined for those values of frequency equal to any of the 
natural frequencies of the component elements* This fact . 
which ^has been illustrated for an axially loaded beam is 
equally valid for torsional or Jlexural modes of vibration. 
)Furthermor^, it should be emphasized that any of these critical 
values (the natural frequencies for isolated members) actually 
may or may not be a natural frequency for the structure as a 
whole. 

The nature of the difficulty in finding the natural fre- 
quencies by equating the determinant of the dynamic stiffness 
'matrix to zero lies in the mathematical condition (equation (9) 
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for axial vibration) required to obtain the dynamic stiffness 
matrix of ihe element. FailOre to give due mathematical regard 
to conditions of the nature of equation (9) may produce either 
a spurious solution or a paradoxical situation such as the one 
described in this article. 

Nov/ make the determinant equal to infinity ! 

Another situation in which lack of a correct and careful 
mathematical treatment leads at best to embarrassing results 
is given by the following case described in detail elsewhere^. 

This case deals with the dynamic analysis of continuous 
beams presented in various texts on structural dynamics 
(Timoshenko (3), Biggs (5), F^rtis (6)). The analysis consists 
of establishing the compatibility conditions at the intermediate 
supports of the beam. These conditions lead to an equation 
relating the redundant moments at three consecutive supports, 
of the same form as the well-known Equatioi) of Three Moments 
in statics. Using this procedure, one equation may be written 
for each internal support of the continuous beam; the result 
is^ a system of linear homogeneous equations. In order for any 
free vibration of the beam to be possible, the determinant of ( j 
the coefficients of the system of equations must be equal to 
zero. _ ' 

To illustrate this method, the particular case of a con- 
tinuous beam of two equal spans invariably is presented in 

* 
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texts' on structural dynamics. For this case there is one 
redundant moment and only one equation. The corresponding 
coefficient of this equation is set equal to zero, resulting 
in the equation 

2{cothA^L - cotx ^L) = 0* (11) 

where 



\ = . (12) 

In this equation, L is the length of one span of the beam, 

is the natural frequency corresponding to the n-mode, m the 

mass of the beam per unit of length and EI its stiffness. 

The roots of equation (11) are then determined numerically, or 

graphically by plotting y^ = cothx^ and yg = cotX^L as functions 

of the argument x^L. The first few of these roots are thus 

found to be t 



X ^L'= 3.92, 7.06, 10.2,^ 
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Up to this point the procedure is a straightforv/ard method 
for finding the eigenvalues o£ the sys"tem and the results 
appear to be satisfactory. However, as it is correctly realized 
by the authors of the texts on this subject, a whole series of 
eigenvalues are not included among the roots listed above. 
According to these authors, the missing values are then deter- 
mined by letting, th'e expression on the left side of eqpaation (11) 
be equal to infinity, that is 

cothx - cot x^L = 1 «> • > ' • (13) 

11 .... ■ 
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which then gives the missing roots as 

X^L =^ TT, 27r, Stt, 

The natural frequencies may then.be found by substituting 
into equation (12) values of X^^ and solving for oj Although 
the final results^ give the correct numerical values for the 
na''Sural frequencies of the two spans beam, it leaves the readpr 
in a quandary with no explanation forthcoming for the formula- 
tion of equation (13 )» in which the left side of equation (11) 
was set equal to - oo, ' 

The elucidation of the answer to this question requires 
* a careful examination of the mathematical development to the 
classic Equation of Three moments for vibration of continuous 
beams. In this derivation it is necessary to divide by the 

factor sin sL, v;hich for this purTDOse should be assumed as 
/ in 

not equal to zero, thus, 

sinXj^L 0 ' (14) 

ConsequEhtiy , this assumption precludes, as possible 
eigenvalues, the roots of sinx^L = 0, which for the particular 
case of a continuous beam of two equal spans, are al-so eiigen- 
values of the p3?oblem. The troublesome, unexplained use of 
equation (13) can be eliminated ^if the analysis is carried out 
mathematically for the sr>ecial case of the beam of two equal 
spans. For this case, the reduction and expansion of the • 
determinantat equation is given by 

12 ' ■ 
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'■. ■ ' ^ 

sin-^ X^L sinh X^L (cos'^X^L sinh X^L - sin x^L cosh X^^L) - 0 (15) 

\ 

from which the characterictic roots are found by setting the 

factors' in equation (15) equal to zero, namely, 

. sinh =0 . . " 

sin XL = X) ' . . , . (17) 

n *> 

cos X^L sinh X - cosh x^L sin x = 0 , (l8.) 

Equation (l6) gives ^nly the trivial solution x^L = 0, while 

equation (l?) gives correctly the series of eigenvalues 

7r,27r,37T.. and equation (l8) which may be written as 

equation (11 ) for siyj^r^L ^ 0 and sinhx 0 gives, as before, 

the series of eigenvalues 3.92, 7.06, 10.2 

» 

• •« • * 

Conclusions ^ 

Three cases have been presented to illustrate that the 
lack of sound mathematical treatment of an engineering problem 
could result in an unacceptable or erroneous analysis* Many 
other examples from the ^ classroom or professional practice may 
be added to fthe cases pi^esented. I'^/ithout taking an extreme 
position, I am inclined to believe that engineering students 
should not be exposed exclusively to applied or "practical 
mathematics" but also to rigorous presentations by competent 
professors of mathematics. In this way the student would 
develop an appreciation and respect for mathematical formulations, 
conditions find proofs. 
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